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Abstract — This paper focuses on the 1-to-A broadcast packet 
erasure channel (PEC), which is a generalization of the broadcast 
binary erasure channel from the binary symbol to that of 
arbitrary finite fields GF(g) with sufficiently large q. We consider 
the setting in which the source node has instant feedback of 
the channel outputs of the K receivers after each transmission. 
Such a setting directly models network coded packet transmission 
in the downlink direction with integrated feedback mechanisms 
(such as Automatic Repeat reQuest (ARQ)). 

The main results of this paper are: (i) The capacity region 
for general l-to-3 broadcast PECs, and (ii) The capacity region 
for two classes of 1-to-A' broadcast PECs: the symmetric PECs, 
and the spatially independent PECs with one-sided fairness con- 
straints. This paper also develops (iii) A pair of outer and inner 
bounds of the capacity region for arbitrary 1-to-A' broadcast 
PECs, which can be evaluated by any linear programming solver. 
For most practical scenarios, the outer and inner bounds meet 
and thus jointly characterize the capacity. 

Index Terms — Packet erasure channels, broadcast capacity, 
channel output feedback, network code alignment. 

I. Introduction 

In the last decade, the new network coding concept has 
emerged (Toll, which focuses on achieving the capacity of a 
communication network. More explicitly, the network-coding- 
based approaches generally model each hop of a packet- 
based communication network by a packet erasure channel 
(PEC) instead of the classic Gaussian channel [3). Such simple 
abstraction allows us to explore the information-theoretic 
capacity of a much larger network with mathematical rigor 
and also sheds new insights on the network effects of a 
communication system. One such example is the broadcast 
channel capacity with message side information. Unlike the 
existing Gaussian Broadcast Channel (GBC) results that are 
limited to the simplest 2-user scenario 1171 . the capacity region 
for 1-to-A broadcast PECs with message side information has 
been derived for K = 3 and tightly bounded for general K 
values lfT31 . 1161 Fl In addition to providing new insights on 
network communications, this simple PEC-based abstraction 
in network coding also accelerates the transition from theory 
to practice. Many of the capacity-achieving network codes |6] 
have since been implemented for either the wireline O or the 
wireless multi-hop networks Q, GO- 

Motivated by recent wireless network coding protocols, this 
paper studies the memory less 1-to-A' broadcast PEC with 

'The results of l-to-K broadcast PECs with message side information 1151 . 
1161 is related to the capacity of the wireless "XOR-in-the-air" scheme (7). 



Channel Output Feedback (COF). Namely, a single source 
node sends out a stream of packets wirelessly, which carries 
information of K independent downlink data sessions, one 
for each receiver dk, k = 1, ■•■ ,K, respectively. After 
packet transmission through the broadcast PEC, each dk then 
informs the source its own channel output by sending back 
the ACKnowledgement (ACK) packets after each time slot. 
derives the capacity region of the memoryless l-to-2 broadcast 
PEC with COF. The results show that COF strictly improves 
the capacity of the memoryless l-to-2 broadcast PEC, a 
mirroring result to the achievability results of GBCs with COF 
fl2l . Other than increasing the achievable throughput, COF 
can also be used for queue and delay management ifTTI . |[T4ll 
and for rate-control in a wireless network coded system (HI- 

The main contribution of this work includes: (i) The capac- 
ity region for general l-to-3 broadcast PECs with COF; (ii) 
The capacity region for two classes of 1-to-A broadcast PECs 
with COF: the symmetric PECs, and the spatially independent 
PECs with one-sided fairness constraints; and (iii) A pair of 
outer and inner bounds of the capacity region for general 1-to- 
K broadcast PECs with COF, which can be evaluated by any 
linear programming solver. Extensive numerical experiments 
show that the outer and inner bounds meet for almost all 
practical scenarios and thus effectively bracket the capacity. 

The capacity outer bound in this paper is derived by 
generalizing the degraded channel argument first proposed in 
1121 . For the achievability part of (i), (ii), and (iii), we devise 
a new class of inter-session network coded schemes, termed 
the packet evolution method. The packet evolution method is 
based on a novel concept of network code alignment, which 
is the PEC-counterpart of the interference alignment method 
originally proposed for Gaussian interference channels Q], Q|. 

This paper is organized as follows. Section [n] contains 
the basic setting and the detailed comparison to the existing 
results. Section ITTT1 describes the main theorems of this paper. 
Section [IV] provides detailed description of the packet evolu- 
tion scheme and the corresponding intuitions. Section HVl also 
includes brief sketches on how to use the packet evolution 
method to prove the achivability results. (Most proofs of 
this paper are omitted due to the limit of space.) Numerical 
evaluation is included in Section [V] Section [VTl concludes this 
paper. 



II. Problem Setting & Existing Results 
A. The Memoryless 1-to-K Broadcast Packet Erasure Channel 

For any positive integer A, we use [A] = {1, 2, • • • , A} to 
denote the set of integers from 1 to A", and use 2^1 to denote 
the collection of all subsets of [A]. 

Consider a 1-to-A" broadcast PEC from source s to K 
destinations dk, k G [AT]. For each channel usage, the 1-to- 
K broadcast PEC takes an input symbol Y G GF(g) from 

s and outputs a A'-dimensional vector Z = {Z\, ■ ■ ■ , Zk) £ 
({y}U{*}) K , where the fc-fh coordinate Zk being "*" denotes 
that the transmitted symbol Y does not reach the fc-th receiver 
dk (thus being erased). There is no other type of noise, i.e., the 
individual output is either equal to the input Y or an erasure 
"*." The success probabilities of a 1-to-A" PEC are described 
by 2 K non-negative parameters: P s r^vrg for all S G 2^ such 
that Ese2iK] P S [k]\s = 1 and for a11 V e GF (<?), 

Prob ({ke[K]:Z k = y} = S\Y = y)=p sWK g. 

That is, P gTyjrg denotes the probability that the transmitted 
symbol Y is received by and only by the receivers {dk : k G 
5}. For all 5 G 2^, we also define 
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That is, pug is the probability that at least one of the receiver 
dk in S successfully receives Y. We sometimes use p k as 
shorthand for Pu{k}> which is the marginal probability that 
the fc-th receiver dk receives Y successfully. 

We assume that the broadcast PEC is memoryless and time- 
invariant, and use Y(t) and Z(t) to denote the input and 
output for the t-th time slot. Note that this setting allows the 
success events among different receivers to be dependent, also 
defined as spatial dependence. For example, when two logical 
receivers dk 1 and dk 2 are situated in the same physical node, 
we simply set the Ps \k\\S P arameters to allow perfect corre- 
lation between the success events of dk 1 and dk 2 . Throughout 
this paper, we consider memoryless 1-to-A broadcast PECs 
that may or may not be spatially dependent. 

B. Broadcast PEC Capacity with Channel Output Feedback 

We consider the following broadcast scenario from s to 
{dk ■ Vfc G [AT]}. Assume slotted transmission. Source s is 
allowed to use the 1-to-A" PEC exactly n times and would 
like to carry information for A" independent downlink data 
sessions, one for each dk, respectively. For each fc G [A], the 
fc-th session (from s to dk) contains nRk information symbols 
X fe = {X k ,j G GF(g) : Vj G [nR k ]}, where R k is the data 
rate for the (s,dk) session. All the information symbols Xkj 
for all fc G [A] and j G [nRk] are independently and uniformly 
distributed in GF(q). 

We consider the setting with instant channel output feedback 
(COF). That is, for the t-th time slot, s sends out a symbol 

Y(t) = f t ({X k : Vfc G [A]}, {Z(r) : r G [t - 1]}) , 



which is a function / t (-) based on the information symbols 
{Xk.j} and the COF {Z(r) : r G [t - 1]} of the previous 
transmissions. In the end of the n-th time slot, each dk outputs 
the decoded symbols 

±k = {X^ : Vj G [nR k ]} = 9 k{{Z k {t) : Vt G [n]}), 

where gk{-) is the decoding function of dk based on the cor- 
responding observation Zk (t) for all t G [n] . Note that we as- 
sume that the PEC channel parameters {ffg pq\s : VS G 2^1 1 
are available at s before transmission. 

We now define the achievability of a 1-to-A PEC with COF. 

Definition 1: A rate vector (Ry, ■ ■ ■ ,Rk) is achievable if 
for any e > 0, there exist sufficiently large n and sufficiently 
large underlying finite field GF(g) such that 



Vfc G [A], Prob(X fc ^X fc J < e. 

Definition 2: The capacity region of a 1-to-A" PEC with 
COF is the closure of all achievable rate vectors. 

C. Existing Results 

Theorem 1 (Theorem 3 in H5]/): The capacity region 
(Ri,R 2 ) of a l-to-2 PEC with COF is described by 
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One scheme that achieves the above capacity region in (Q]) 
is the 2-phase approach in |5|. That is, for any (Ai,A2) in 
the interior of ((T), perform the following coding operations. 

In Phase 1, s sends out uncoded information packets Xxj 1 
and X2j 2 for all j\ G [nRy] and j2 £ [nR-2\ until each packet 
is received by at least one receiver. Those Xij 1 packets that 
are received by d\ have already reached their intended receiver 
and thus will not be retransmitted in the second phase. Those 
Xij 1 packets that are received by d^ but not by d\ need to 
be retransmitted in the second phase, and are thus stored in 
a separate queue Q 1 . 2 j- Symmetrically, the X2j 2 packets that 
are received by d% but not by d 2 need to be retransmitted, and 
are stored in another queue Q 2 -i2- Since those "overheard" 
packets in queues Q 1 . 2 j and Q 2 12 perfect candidates for 
intersession network coding Q, they can be linearly mixed 
together in Phase 2. Each single coded packet in Phase 2 can 
now serve both d\ and d 2 simultaneously. The intersession 
network coding gain in Phase 2 allows us to achieve the 
capacity region in fl}. Based on the same logic, J9) derives 
an achievability region for 1-to-A" broadcast PECs with COF 
under a perfectly symmetric setting. iTOIl implements such 2- 
phase approach while taking into account of various practical 
considerations, such as time-out and network synchronization. 

D. The Suboptimality of The 2-Phase Approach 

The above 2-phase approach does not achieve the capacity 
for the cases in which K > 2. To illustrate this point, consider 
the example in Fig. [T] 

In Fig. Q] source s would like to serve three receivers d\ to 
g?3. Each (s,dk) session contains a single information packet 
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Fig. 1 . Example of the suboptimality of the 2-phase approach. 



Xk, and the goal is to convey each Xk to the intended dk for 
all A; = 1,2, 3. Suppose the 2-phase approach in Section Hl-CI 
is used. During Phase 1, each packet is sent repeatedly until 
it is received by at least one receiver, which either conveys 
the packet to the intended receiver or creates an overheard 
packet that can be used in Phase 2. Suppose after Phase 1, d\ 
has received X2 and X%, d2 has received Xi and X3, and d 3 
has not received any packet (Fig. [TJ). Since each packet has 
reached at least one receiver, source s moves to Phase 2. 

Suppose s sends out a coded packet [Xi + X%] in Phase 2. 
Such coded packet can serve both d\ and g?2- That is, d\ 
(resp. 1^2) can decode X\ (resp. X2) by subtracting X2 (resp. 
Xi) from [X\ +X2]. Nonetheless, since the broadcast PEC is 
random, the packet [Xi +X2] may or may not reach di or c?2- 
Suppose that due to random channel realization, [X% + X2] 
reaches only cfe, see Fig.Q] The remaining question is what s 
should send for the next time slot. 

The existing 2-phase approach: We first note that since 
c?3 received neither Xi nor X2 in the past, the newly received 
[Xi + X2] cannot be used by d$ to decode any information 
packet. In the existing results Q, (9), d 3 thus discards the 
overheard [Xi + X2], and s would continue sending [Xi + 
X2] for the next time slot in order to capitalize this coding 
opportunity created in Phase 1. 

The optimal decision: It turns out that the broadcast system 
can actually benefit from the fact that d 3 overhears the coded 
packet [Xi + X2] even though neither X\ nor X2 can be 
decoded by cfo. More explicitly, instead of sending \X\ +X2], 
s should send a new packet [X% + X2 + X3] that mixes all 
three sessions together. With the new [X\ + X2 + X3] (plus 
the previous overhearing patterns in Fig. [TJ, d\ can decode 
Xi by subtracting both X 2 and X 3 from [Xi + X 2 + X 3 }. 
d 2 can decode X 2 by subtracting both X\ and X 3 from 
[Xi + X 2 + X 3 ]. For d 3 , even though da does not know the 
values of X\ and X2, d-$ can still use the previously overheard 
\X\ +X2] packet to subtract the interference [X\ + X2) from 
[X\ + X2 + X3} and decode its desired packet X3. As a 
result, the new coded packet [X\ + X2 + X3] serves d\, d2, 
and ds, simultaneously. This new coding decision thus strictly 
outperforms the existing 2-phase approach. 

Two critical observations can be made for this example. 
First of all, when d 3 overhears a coded \X\ + X2] packet, 
even though 1^3 can decode neither X\ nor X2, such new 
side information can still be used for future decoding. More 
explicitly, as long as s sends packets that are of the form 
ct{Xi + X 2 ) + PX 3 , the "aligned interference" a(X 1 + X 2 ) 



can be completely removed by d 3 without decoding individual 
X\ and X2. This technique is thus termed "code alignment, 7 ' 
which is in parallel with the original interference alignment 
method [fl~). Second of all, in the existing 2-phase approach, 
Phase 1 has the dual roles of sending uncoded packets to their 
intended receivers, and, at the same time, creating new coding 
opportunities (the overheard packets) for Phase 2. It turns out 
that this dual-purpose Phase- 1 operation is indeed optimal. 
The suboptimality of the 2-phase approach for K > 2 is 
actually caused by the Phase-2 operation, in which source s 
only capitalizes the coding opportunities created in Phase 1 but 
does not create any new coding opportunities for subsequent 
packet mixing. One can thus envision that for the cases K > 2, 
an optimal policy should be a multi-phase policy, say an M- 
phase policy, such that for all i G [M — 1] (not only for 
the first phase) the coded packets sent in the i-th phase have 
dual roles of carrying information to their intended receivers 
and simultaneously creating new coding opportunities for the 
subsequent Phases (i + 1) to M. These two observations will 
be the building blocks of our achievability results. 

III. The Main Results 



Section IIII-AI focuses on the capacity results for arbitrary 
broadcast PEC parameters while Section IIII-BI considers two 
special classes of broadcast PECs: the symmetric and the 
spatially independent PECs, respectively. 

A. Capacity Results For General 1-to-K Broadcast PECs 

We define any bijective function tt : [K] \-> [K] as a 
permutation. There are totally K\ distinct permutations. Given 
any permutation n, for all j e [K] we define SJ = {tt(1) : 
V/ G [j]} as the set of the first j elements according to the 
permutation w. We then have the following capacity outer 
bound for any l-to-K broadcast PEC with COF. 

Proposition 1: Any achievable rates {R\,--- ,Rk) must 
satisfy the following K\ inequalities: 



K R 
W , £fMzl<i. 
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Sketch of the proof: For any given tt, construct a new 
broadcast channel from the original one by adding (K — 1) 
information pipes connecting all the receivers d\ to dn ■ More 
explicitly, for all j G [K — 1], create an auxiliary pipe from 
dn(j) to d v (j+iy With the new auxiliary pipes, the success 
probability of d„(j) increases from p n u) to pus™ f° r all j G 
[K] since d^U) now knows the transmitted symbol Y as long 
as at least one of d^m, V/ G [j], receives Y successfully. Note 
that the new broadcast PEC is physically degraded. By the 
same arguments as in Q, |fl~2], (0 describes the capacity of a 
physically degraded PEC with COF, which thus outer bounds 
the capacity of the original PEC with COF. 

■ 

For the following, we provide the capacity results for 
general l-to-3 broadcast PECs. 



Proposition 2: For any parameter values 
{ps {i 2 3 }\s : V5 e 2^' 2 ' 3 }| of a l-to-3 PEC with COF, the 
capacity outer bound in Proposition Q] is the capacity region. 

To state the capacity inner bound for 1-to-A PECs with 
K > 4, we need to define an additional function: f p (ST), 
which takes an input ST of two disjoint sets S, T G 2^1. 
More explicitly, f p (ST) is the probability that the transmitted 
packet Y is received by all those di with i £ S but not received 
by any dj with j 6 T. That is, 



f P (ST) ± 
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We also say that a ifn'cr foto/ ordering " on 2^1 is 
cardinality-compatible if 

V5 , 1 ,5 2 e2^, [S x | < \S 2 \ => Si -< S 2 . 

For example, for A" = 3, the following strict total ordering 

-< {2} -< {1} -< {3} -< {1,2} -< {1,3} -< {2,3} -< {1,2,3} 

is cardinality-compatible. 

Proposition 3: Fix any cardinality-compatible, strict total 
ordering -<. For any l-to-K PEC with COF, a rate vector 
(Ri, • • ■ , Ar-) can be achieved by a linear network code if 
there exist 2 K non-negative x variables, indexed by S G 2^: 



|x s >0:VS'G2[ A T}, 



(3) 



and A"3 X 1 non-negative w variables, indexed by (fc; S — > T) 
satisfying T C 5 C ([A]\fc): 

{w* ; s->t > : Vfc e [A], VS,T 6 2^ 

satisfying T C S C ([A"]\fc)} , (4) 
such that jointly the following linear inequalities^] are satisfied: 

(5) 



E x * < 1 

VS:S<E2M 

VT e 2^, Vfc e r, 

X7 1 > 



E W k ;S^(T\k) (6) 

VS:(T\fc)CSC([K]\fc) 

Vfc G [A"], w fc; 0^0 -pu[Jf] > Rk (7) 
Vfce [A],VSC ([A]\fc),S/0, 

^2 Wk-S-^Tt Pu([K]\S) > 
,VTi:7iCS / 

E 



"Mi-^ • f P (S\T X )([A]\5) 



VSi , Ti : such that 
TiCSiC ([if]\fc), 
Tj C S, 5 £ Si 



(8) 



2 There are totally (1 + K2 K_1 +ii"3^ _1 ) inequalities. More explicitly, 
|5J describes one inequality. There are K2 K ~ 1 inequalities having the form 
of (6). There are totally K3 K_1 inequalities having the form of one of 0, 
{sj, and (9). For comparison, the outer bound in Proposition ff] actually has 
more inequalities asymptotically (K\ of them) than those in Proposition fj] 



Vfc G [A], S, T e 2 [K] satisfying T C S C ([K]\k),T ^ S, 
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VTi C S : 
(TiU{fc}) -< (Tu{fe}) 
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VSi : Si -< S, 
T C S x C ([K]\*0 



E 



VSi , Ti : such that 

Ti C Si C ([if]\fc), 
(Ti U {k}) X(TU {fe}), 
Ti C S, S £ Si 



^;Si^Ti •/ P ((S\T 1 )([A]\5)) 



(9) 



Remark: For some general classes of PEC parameters, one 
can prove that the inner bound of Proposition [3] is indeed the 
capacity region for arbitrary A" > 4 values. Two such classes 
are discussed in the next subsection. 

B. Capacity Results For Two Classes of l-to-K PECs 

We first focus on symmetric broadcast PECs. 
Definition 3: A 1-to-A broadcast PEC is symmetric if the 
channel parameters {pg T^]\g ; ^S G 2^1 j satisfy 

VS lf S 2 G 2^1 with \Si\ = |S 2 |, P SlW ]^ = P S2Wm . 

Proposition 4: For any symmetric 1-to-A" broadcast PEC 
with COF, the capacity outer bound in Proposition Q] is indeed 
the corresponding capacity region. 

In addition to perfect channel symmetry, another practical 
setting is to allow channel asymmetry while assuming spatial 
independence between different destinations di. 

Definition 4: A 1-to-A broadcast PEC is spatially indepen- 
dent if the channel parameters |pg Tg]\g : ^S G 2^1 1 satisfy 



VS G 2^ K \ p 



S[K]\S 



\kes 



n (i-ffc) 

y ke[K]\S 



where pk is the marginal success probability of destination dk- 

To describe the capacity results for spatially independent 
1-to-A PECs, we need the following additional definition. 

Definition 5: Consider a 1-to-A" broadcast PEC with 
marginal success probabilities p\ to pk- We say a rate vector 
• • • , Rk) is one-sidedly fair if Vi ^= j satisfying pi < pj, 
we have Ri(l — pi) > Rj(l —Pj)- We use A os f to denote the 
collection of all one-sidedly fair rate vectors. 

The one-sided fairness contains many practical scenar- 
ios of interest. For example, the perfectly fair rate vector 
(R, R, ■ ■ ■ , R) by definition is also one-sidedly fair. Another 
example is when min(pi, • • • ,px) > |> a proportionally fair 
rate vector (piA,p 2 A, • ■ • ,pkR) is also one-sidedly fair. 

For the following, we provide the capacity of spatially 
independent 1-to-A" PECs with COF under the condition of 
one-sided fairness. 



Proposition 5: Suppose the 1-to-A' PEC of interest is spa- 
tially independent and the marginal success probabilities sat- 
isfy < pi < P2 < • • • < Pk, which can be achieved by 
relabeling. Any (iii, ■ ■ ■ , Rk) £ A os f is in the capacity region 
if and only if (Ri, ■ ■ ■ , Rk) G A os f satisfies 
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R k 



K 

y 



< i. 



(10) 



Namely, Proposition Q] is indeed the capacity region when 
focusing on the one-sidedly fair rate region A 0S f. 

IV. The Packet Evolution Schemes 

We now describe a new class of coding schemes, termed 
the packet evolution (PE) scheme, which is the building block 
of the capacity / achievability results in Section [HI] 

A. Description Of The Packet Evolution Scheme 

Recall that each (s,dk) session has nRk information pack- 
ets Xk.i to Xk.nR k - We associate each of the X^fcLi n Rk 
information packets with an intersession coding vector v and a 
set S C [A]. An intersession coding vector is a ^X^tLi n Rkj - 
dimensional row vector with each coordinate being a scalar in 
GF(g). Before the start of the broadcast, for any k £ [A'] 
and j £ [nRk] we initialize the corresponding vector v of 
X k ,j in a way that the only nonzero coordinate of v is the 
coordinate corresponding to X k j and all other coordinates 
are zero. Without loss of generality, we set the value of the 
only non-zero coordinate to one. That is, initially the coding 
vectors v are set to the elementary basis vectors. 

For any k £ [A] and j £ [nRk] the set S of X k j is 
initialized to 0. We call S the overhearing set of the packet 
X k j. We use v(X kj j) and S(Xk.j) to denote the intersession 
coding vector and the overhearing set of a given X kt j, 

Throughout the n broadcast time slots, s constantly updates 
S(Xk.j) and v(X k j) according to the COF. The main struc- 
ture of a packet evolution scheme can now be described as 
follows. 

§ The Packet Evolution Scheme 
1: Source s maintains a flag f change- Initially, set flange <— E 

2: for t = 1, • • ■ ,n, do 

3: In the beginning of the i-th time slot, do Lines l4l to [TOl 

4: if fchange = 1 then 

5: Choose a non-empty subset T C [A]. 

6: Run a subroutine Packet Selection, which takes 
T as input and outputs a collection of |T| packets 
{Xk,j k ■ Vfc £ T}, termed the target packets. The 
output {X kJk } must satisfy (S(X kJk ) U {k}) D T 
for all k £ T. 

7: Generate k uniformly random coefficients c k £ 
GF(g) and construct an intersession coding vector 

8: Set fchange ^ 0. 



end if 

Sends out a linearly intersession coded packet according 
to the coding vector v tx . That is, we send 



Y u — Vtx ■ (-^1,1 j 



,X 



K.nRj 



where (Xi t i, ■ ■ ■ ,Xfc, n R K ) T i s a column vector con- 
sisting of all information symbols. 

11: In the end of the t-th time slot, use a subroutine 
Update to revise the v(X kjk ) and S(X k , Jk ) values 
of all target packets X k ,j k based on the COF. 

12: if the S(X k j k ) value changes for at least one target 
packet Xk,j h after the Update then 

13: Set fchange <~ E 

14: end if 
15: end for 



In summary, a group of target packets {X k ,j k } are selected 
according to the choice of the subset T. The corresponding 
vectors {v(X k j k )} are used to construct a coding vector v tx . 
The same coded packet Y tx , corresponding to v tx , is then sent 
repeatedly for many time slots until one of the target packets 
X k ,j k evolves (when the corresponding S(X k ,j k ) changes). 
Then a new subset T is chosen and the process is repeated 
until we use up all n time slots. Three subroutines are used 
as the building blocks of a packet evolution method: (i) How 
to choose the non-empty T C [A']; (ii) For each k £ [A], 
how to select a single target packets X k ,j k among all X k j 
satisfying (S(X k ,j) U {k}) D T; and (iii) How to update the 
coding vectors v(X k j k ) and the overhearing sets S(X k .j k ). 
We first describe the detailed update rule of (iii). 

§ Update of S(X kijk ) and v(X kJk ) 
1: Input: The T and v tx used for transmission in the current 
time slot; And S m , the set of destinations d$ that receive 
the transmitted coded packet in the current time slot. 



for all k £ T do 

if S a £S(X k>jh ) then 

Set S{X k . 3k ) <- (Tn S(X kijk )) U S B 
Set v(X k j k ) <— vtx. 
end if 
end for 



An Illustrative Example Of The PE Scheme: 

Let us revisit the optimal coding scheme of the example 
in Fig. Q] of Section III-DI After initialization, the three infor- 
mation packets X\ to A3 have the corresponding v and S: 

v(JTi) = (1,0,0), v(A 2 ) = (0,1,0), and v(A 3 ) = (0,0,1), 
and S(X{) = S(X 2 ) = S(X 3 ) = 0. We use the following 
table for summary. 



Xi. (1,O,O),0 A 2 : (O,1,O),0 X 3 : (0,0,1), 



Consider a duration of 5 time slots. 

Slot 1: Suppose that s chooses T = {1}. Since (0U{1}) D 
T, Packet Selection outputs X x . The coding vector v tx 
is thus a scaled version of v(Ai) = (1,0,0). Without loss of 



generality, we choose v tx = (1,0,0). Based on v tx , s transmits 
a packet lXi + 0X 2 + 0X 3 = X\. Suppose [Xi] is received 
by d 2 , i.e., S a = {2}. Then during UPDATE, S a = {2} £ 
5(Xi) = 0. Update thus sets S(X{) = {2} and v(X x ) = 
v te = (1)0,0). The packet summary becomes 



X x : (1,0,0),{2} X 2 : (O,1,O),0 X 3 : (0,0,1) 



Slot 2: Suppose that s chooses T = {2}. Since (0U{2}) D 
T, Packet Selection outputs X 2 . The coding vector v tx is 
thus a scaled version of v(X 2 ) = (0,1,0). Without loss of 
generality, we choose v tx = (0, 1,0) and accordingly [X 2 ] is 
sent. Suppose [X 2 ] is received by d\, i.e., S IX = {1}. Since 
Srx ^ S(X 2 ), after UPDATE the packet summary becomes 



Xr. (1,0,0),{2} X 2 : (0,1,0),{1} X 3 : (0,0,1), 



Slot 3: Suppose that s chooses T = {3} and PACKET 
Selection outputs X 3 . v tx is thus a scaled version of 
v(X 3 ) = (0, 0, 1), and we choose v tx = (0, 0, 1). Accordingly 
[X3] is sent. Suppose [X3] is received by {di,d 2 }, i.e., 
Si% = {1, 2}. Then after UPDATE, the summary becomes 



X x : (1,0,0),{2} X 2 : (0,1,0),{1} X 3 : (0,0,1),{1, 2} 



Slot 4: Suppose that s chooses T = {1, 2}. Since (S(Xi)U 
{1}) 2 T and (S{X 2 ) U {2}) D T, Packet Selection 
outputs {Xi, X 2 }. Vtx is thus a linear combination of v(X\) = 
(1, 0, 0) and v(X 2 ) = (0, 1, 0). Without loss of generality, we 
choose Vtx = (1,1,0) and accordingly [X\ + X 2 ] is sent. 
Suppose [X\ + X 2 ] is received by d 3 , i.e., S rx = {3}. Then 
during Update, for X u S a = {3} £ S(Xi) = {2}. Update 
thus sets = {2,3} and v(X{) = v tx = (1,1,0). For 

X 2 , S a = {3} £ S(X 2 ) = {1}. Update thus sets S(X 2 ) = 
{1,3} and v(X 2 ) = v tx = (1, 1,0). The summary becomes 



Xf. (1,1,0),{2,3} 


X 2 : (1,1,0),{1,3} 


X 3 : (0,0,1),{1,2} 





Slot 5: Suppose that s chooses T = {1, 2, 3}. By Line|6]of 
The Packet Evolution Scheme, the subroutine Packet 
Selection outputs {Xi,X 2 ,X 3 }. v tx is thus a linear com- 
bination of v(Xi) = (1,1,0), v(X 2 ) = (1,1,0), and 
v(X 3 ) = (0, 0, 1), which is of the form a(X 1 + X 2 ) + j3X 3 . 
Note that the packet evolution scheme automatically achieves 
code alignment, which is the key component of the optimal 
coding policy in Section III-DI Without loss of generality, we 
choose a = /3 = 1 and v tx = (1, 1, 1). Y tx = [Xi + X 2 + X 3 ] 
is sent accordingly. Suppose \X\ + X 2 + X 3 ] is received by 
{di,d 2 ,d 3 }, i.e., S IX = {1,2,3}. Then after UPDATE, the 
summary of the packets becomes 



X x : (1,1,1),{1,2,3} 


X 2 : (1,1,1),{1,2,3} 


X 3 : (1,1,1),{1,2,3} 





From the above step-by-step illustration, we see that the 
optimal coding policy in Section III-DI is a special case of a 
packet evolution scheme. 

B. Properties of A Packet Evolution Scheme 

We term the packet evolution (PE) scheme in Section lTV-Al a 
generic PE method since it does not depend on how to choose 



T and the target packets Xk,j k and only requires the output of 
Packet Selection satisfying (S(X k . jk ) U {k}) D T,Vfc e 
T. In this subsection, we state some key properties of any 
generic PE scheme. The intuition of the PE scheme is based 
on these key properties and will be discussed in Section HV-CI 

We first define the following notation for any linear network 
codes. (Note that the PE scheme is a linear network code.) 

Definition 6: Consider any linear network code. For any 
destination dk, each of the received packet Zk(t) can be 
represented by a vector Wfc(t), which is a ^X^tLi n Rkj - 
dimensional vector containing the coefficients used to generate 
Z fc (t). That is, Z k {t) = w fe (t) • {X ltU --- ,X K<nRK ) T . If 
Zk{t) is an erasure, we simply set Wfc(i) to be an all-zero 
vector. The knowledge space of destination dk in the end of 
time t is denoted by Qz,fc(t), which is the linear span of Wfc(r), 
t <t. That is, fi Zifc (t) = span(w fe (r) : Vr 6 [t]). 

Definition 7: For any non-coded information packet 
Xk,j, the corresponding intersession coding vector is a 
(^2^ =1 nRkj -dimensional vector with a single one in the 
corresponding coordinate and all other coordinates being zero. 
We use 8kj to denote such a delta vector. The message space 
of dk is then defined as fiju.fc = span(<5fcj : Vj £ [nRk]). 

The above definitions imply the following straightforward 
lemma: 

Lemma 1: In the end of time t, destination dk is able to 
decode all the desired information packets Xkj, Vj G [nRk], 
if and only if Q M ,fe Q foz,k(t). 

We now define "non-interfering vectors" from the perspec- 
tive of a destination dk- 

Definition 8: In the end of time t (or in the beginning of 
time (t + 1)), a vector v (and thus the corresponding coded 
packet) is "non-interfering" from the perspective of dk if 

v G span(n z .k{t),n M ,k)- 



By definition, any non-interfering vector v can always 
be expressed as the sum of two vectors v' and w, where 
v' G Q>M,k is a linear combination of all information vectors 
for dk and w G £lz,k(t) is a linear combination of all the 
packets received by dk- If v' = 0, then v = w is a transparent 
packet from dk's perspective since dk can compute the value 
of w- (X\ t i, • • • , Xk ,nR. K ) T from its current knowledge space 
rtz,k(t)- If v ' 7^ 0, then v = v' + w can be viewed as 
a pure information packet v' G Qm,*; after subtracting the 
unwanted w vector. In either case, v is not interfering with 
the transmission of the (s,dk) session, which gives the name 
of "non-interfering vectors." 

The following Lemmas|2]and[3]discuss the time dynamics of 
the PE scheme. To distinguish different time instants, we add 
a time subscript and use St-i(Xk j k ) and St(Xk.j k ) to denote 
the overhearing set of Xk,j k in the end of time (t—1) and t, 
respectively. Similarly, v t -i(Xkj k ) and v t (Xkj k ) denote the 
coding vectors in the end of time (t — 1) and t, respectively. 

Lemma 2: In the end of the t-th time slot, consider any 
Xkj out of all the information packets X 1 1 to Xk, u r k - Its 



assigned vector v t (X k j) is non-interfering from the perspec- 
tive of di for all i e (S t {X kJ ) U {k}). 

To illustrate Lemma|2] consider our 5-time-slot example. In 
the end of Slot 4, we have v(Xi) = (1, 1,0) and U 
{1} = {1,2,3}. It can be easily verified by definition that 
v(Xl) = (1, 1,0) is non-interfering from the perspectives of 
di, d,2, and d%, respectively. 

Lemma 3: In the end of the i-th time slot, we use fln ik (t) 
to denote the remaining space of the PE scheme: 

fifl,k(t) = 

span(v t (X fcJ ) : Vj G [nR k ] satisfying k £ S t (X kij )). 

For any n and any e > 0, there exists a sufficiently large 
finite field GF(g) such that for all k G [K] and t G [n], 

Prob(span(Oz jfe (t),Qi ii fc(i)) = span(fiz jfe (i), Q M ,fc)) 
> 1 - e. 

Intuitively, Lemma [3] says that if in the end of time 
t we directly transmit all the remaining coded packets 
{v t (X ktj ) : Vj G [nR k } 7 k £ S t (X k j)} from s to d k through a 
noise-free information pipe, then with high probability, d k can 
successfully decode all the desired information packets X k} i 
to Xk, n R k ( see Lemma [TJ by the knowledge space flz,k(t) 
and the new information of the remaining space H,R, k (t). 

C. The Intuitions Of The Packet Evolution Scheme 

Lemmas |2] and [3] are the key properties of a PE scheme. In 
this subsection, we discuss the corresponding intuitions. 

Receiving the information packet X k y. Each informa- 
tion packet keeps a coding vector \r(X k j). Whenever we 
would like to communicate X k j to destination d k , instead 
of sending a non-coded packet X k j directly, the PE scheme 
sends an intersession coded packet according to the coding 
vector v(X k j). Lemma [5] shows that if we send all the 
coded vectors v(X k ,j) that have not been heard by d k (with 
k £ S(X k j)) through a noise-free information pipe, then d k 
can indeed decode all the desired packets X k j with close-to- 
one probability. It also implies, although in an implicit way, 
that once a v(X k ,j ) is heard by d k for some jo (therefore 
k G S(X k ,j )), there is no need to transmit this particular 
v(X k j a ) in the later time slots. Jointly, these two implications 
show that we can indeed use the coded packet v(X k j) as 
a substitute for X k ,j without losing any information. In the 
broadest sense, we can say that d k receives a packet X k j 
if the corresponding v(X k j) successfully arrives d k in some 
time slot t. 

Serving multiple destinations simultaneously by mixing 
non-interfering packets: The above discussion ensures that 
when we would like to send X k j k to d k , we can send a 
coded packet \r(X k j k ) as a substitute. On the other hand, by 
Lemma|2] such \r(X k j k ) is non-interfering from di's perspec- 
tive for all i G (S(X k j k )U{k}). Therefore, instead of sending 
a single packet v(X k j k ), it is beneficial to linearly combine 
the transmission of two packets v(X kt j k ) and v(Xij t ) to- 
gether, as long as I G S(X k j k ) and k G S(Xij l ), Since 



v(X k j k ) is non-interfering from di's perspective, it is as if di 
directly receives v(Xij l ) without any interference. Similarly, 
since vpQjJ is non-interfering from d k 's perspective, it is as 
if d k directly receives v(X k j k ) without any interference. By 
generalizing this idea, the PE scheme first selects a T C [K] 
and then constructs a v tx that can serve all destinations k G T 
simultaneously by mixing the corresponding non-interfering 
vectors. 

Creating new coding opportunities while exploiting the 
existing coding opportunities: As discussed in the example 
of Section III-DI the suboptimality of the existing 2-phase 
approach for K > 3 destinations is due to the fact that it 
fails to create new coding opportunities while exploiting old 
coding opportunities. The PE scheme was designed to solve 
this problem. Let us assume that the PACKET SELECTION in 
Line|6]chooses the X k j such that S(X k j) = T\k. That is, we 
choose the X k j that can be mixed with those (s,di) sessions 
with I G S(X k j) U {k} = T. Then Line g] of the UPDATE 
guarantees that if some other di, i ^ T, overhears the coded 
transmission, we can update S(X k j) with a strictly larger set 
S(X k j) U 5 1X . Therefore, new coding opportunity is created 
since we can now mix more sessions (all di, i G S(X k ,j)) 
together with X k j. Note that the coding vector v(X k j) is 
also updated accordingly. The new v(X k ,j) represents the 
necessary "code alignment" in order to utilize this newly 
created coding opportunity. The (near-) optimality of the PE 
scheme is rooted deeply in the concept of code alignment, 
which aligns the "non-interfering subspaces" through the joint 
use of S(X k j) and v(X k j). 

D. Analysis Of The PE Scheme 

One advantage of a PE scheme is that although different 
packets X k ,j k and Xi ^ with k ^ i may be mixed together, 
the corresponding evolution of X k j k (the changes of S(X k j k ) 
and v(X k j k )) are independent from the evolution of Xij i 
(see Line |2 of the Update). Also by Lemma |2] two different 
packets X k j k and Xij i can share the same time slot without 
interfering each other as long as i G S(X k j k ) and k G 
S(Xi t j i ). As a result, the throughput analysis can be done by 
focusing on the individual sessions separately, and considering 
how many time slots from different sessions can be combined 
together. The achievability results are proven by analyzing 
the throughput of the PE scheme with carefully designed 
mechanisms of choosing the set T and the corresponding target 
packets {X k j k : k G T} of a generic PE scheme. 

V. Numerical Evaluation 

We first notice that both the inner and outer bounds are 
linear programming (LP) problems and can be evaluated 
by any LP solvers. We perform numerical evaluation for 
spatially independent 1-to-K broadcast PECs with COF by 
randomly varying the values of the marginal success proba- 
bilities (pi, - • - ,pk)- Note that although there is no tightness 
guarantee for K > 4 except in the one-sidedly fair rate region, 
in all our numerical experiments with K < 6 (totally 3 x 10 4 of 
them), we have not found any instance of the input parameters 




(a) 



Fig. 2. (a) The sum-rate capacity R* um perf fair in a perfectly fair system versus 
the marginal success probability p of a symmetric, spatially independent 1- 
to-K broadcast PEC, K = 2, 4, 20, and 100. (b) The sum-rate capacities for 
a 6-destination heterogenous channel profiles with the success probabilities 
pi to pe evenly spaced between (p, 1). 

(pi, • • • for which the gap between the outer and inner 

bounds is greater than the numerical precision of the LP solver. 
This shows that Propositions Q] and [3] effectively describe the 
capacity region from the practical perspective. 

To illustrate the network coding gain, we compare the sum- 
rate capacity versus the sum rate achievable by time sharing. 
Fig. |2(a)| considers symmetric, spatially independent PECs 
with marginal success probabilities pi = ■ ■ ■ = pk = P- We 
plot the sum rate capacity R* um perf fair versus p for a perfectly 
fair system R\ = ■ ■ ■ = Rk- As seen in Fig. |2(a)| the network 
coding gains are substantial when we have K > 4 destinations. 
It can also be proven that for any p g (0,1], -R s * um> p er f.fair 
approaches one as K — > oo, which was first observed in J9j- 

We are also interested in the sum rate capacity under asym- 
metric channel profiles (also known as heterogeneous channel 
profiles). Consider asymmetric, spatially independent PECs. 
For each p value, we let the channel parameters p\ to pk be 
equally spaced between (p, 1), i.e., pu = p + (k — l)j^rp 
We then plot the sum rate capacities for different p values. 
(In this experiment, the outer and inner bounds in Section [HI] 
meet for all different p values.) Fig. |2(b)| describes the case for 
K = 6. We plot the curves for perfectly fair (i?i = • • • = Rk) 
and proportionally fair (Rk oc pk) systems, respectively. For 
comparison, we also plot the time-sharing capacity under 
the heterogeneous channel profile. For comparison between 
symmetric (homogeneous) and asymmetric (heterogeneous) 
channel profiles, we plot the sum-rate capacity for symmetric 
channels as well. As can be seen in Fig. |2(b)| network 
coding again provides substantial improvement for all p values. 
However, the gain is not as large as in the case of symmetric 
channels. The results show that for practical implementation, it 
is better to group together all the sessions of similar marginal 
success rates and perform intersession network coding within 
the same group. 

VI. Conclusion 

In this work, we have proposed a new class of intersession 
network coding schemes, termed the packet evolution (PE) 
schemes, for the broadcast PECs with COF. Based on the PE 
schemes, we have derived the capacity region for general 1- 
to-3 broadcast PECs, and a pair of capacity outer and inner 



bounds for general 1-to-A' broadcast PECs, both of which 
can be easily evaluated by any linear programming solver for 
the cases K < 6. It has also been proven that the outer and 
inner bounds meet for two classes of 1-to-A" broadcast PECs: 
the symmetric broadcast PECs, and the spatially independent 
broadcast PECs with the one-sided fairness rate constraints. 
Extensive numerical experiments have shown that the outer 
and inner bounds meet for almost all broadcast PECs encoun- 
tered in practical scenarios. 
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